Spin-like susceptibility of metallic and insulating thin films at low temperature 
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Susceptibility measurements of patterned thin films at sub-K temperatures were carried out using 
a scanning SQUID microscope that can resolve signals corresponding to a few hundred Bohr mag- 
netons. Several metallic and insulating thin films, even oxide-free Au films, show a paramagnetic 
response with a temperature dependence that indicates unpaired spins as the origin. The observed 
response exhibits a measurable out-of-phase component, which implies that these spins will create 
l//-like magnetic noise. The measured spin density is consistent with recent explanations of low 
frequency flux noise in SQUIDs and superconducting qubits in terms of spin fluctuations, and sug- 
gests that such unexpected spins may be even more ubiquitous than already indicated by earlier 
measurements. Our measurements set several constraints on the nature of these spins. 
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The origin of 1//" noise in superconducting quantum 
interference devices (SQUIDs) has been an unresolved 
mystery for more than two decades. Part of this noise 
of typically a few /i^q/VHz at 1 Hz appears to be sur- 
prisingly universal and behaves in every respect like true 
flux noise [l|, Q . Measurements of electron spin dephas- 
ing rates of donors in Si also give evidence for magnetic 
noise originating at or near the surface Recently, 
evidence that a similar noise causes dephasing in super- 
conducting qubits 0, 0, has increased the inter- 
est in this phenomenon. Koch, DiVincenzo and Clarke 
showed that fluctuating electron spins could explain the 
observed magnitude of noise Q . Measurements of a flux 
offset in SQUIDs proportional to 1/T ^ provide direct 
evidence for the presence of spins in superconducting de- 
vices. Even though they are likely related to material 
imperfections such as surface oxides, defect states or con- 
taminations, the nature of these hypothetical spins is in- 
teresting for its own sake, and mitigating their effects is 
essential for several solid-state approaches to quantum 
computation. It is important to understand their origin 
and polarization dynamics, which determine the mag- 
netic noise spectrum. 

Various models assuming different relaxation mecha- 
nisms of unpaired electron spins on defects were recently 
proposed. Koch et al argued that the spin of an electron 
in a charge trap could remain locked until it leaves the 
trap 0. A different model exploring thermally activated, 
nonmagnetic two level systems as cause for spin flips was 
motivated by the argument that only a small fraction of 
all defects have an activation energy low enough to allow 
charge fluctuations [l^. Faoro and loffe explored noise 
from spin diffusion mediated by RKKY coupling via the 
conduction electrons in metallic device elements 



We have measured the magnetic response of metallic 
and insulating thin films from T = 25 mK to 0.6 K, using 
a SQUID susceptometer in a scanning microscope. The 



most prominent result is a surprisingly large paramag- 
netic susceptibility with a 1/T-like temperature depen- 
dence, and a magnitude consistent with a spin- 1/2 den- 
sity of about 4 • 10^^/m^, close to estimates from 1//" 
noise levels in SQUIDs. Furthermore, the response has 
a measurable out-of-phase component, which implies po- 
larization noise through the fluctuation-dissipation theo- 
rem (FDT). Our results thus demonstrate the existence of 
paramagnetic spins, with a density and dynamics suitable 
for producing 1//" noise. Similar susceptibilities seen for 
Ag films and Au films with and without a sticking layer, 
together with previous results inferring the presence of 
spins in superconducting devices, suggest that the spins 
can occur similarly for different materials. 

We focus on results from two samples, which were de- 
signed for other experiments [l2[ and include a range of 
structures with different layer combinations. On sample 
I, rings and wires were e-beam evaporated at a rate of 
about 1.2 nm/s from a 6N purity Au source onto a Si 
substrate with a native oxide [Fig. [11(f), (i)]. First, the 
wires, with widths of 2 and 15 /im, were patterned using 
optical lithography and liftoff. Their thickness was 100 
nm, including a 7 nm Ti adhesion layer. Subsequently, 
the micron-scale, 140 nm thick Au rings, which did not 
include any adhesion layer, were defined using e-beam 
lithography with PMMA [poly(methyl methacrylate)] re- 
sist and liftoff. Some of them were connected to the wider 
wires for heat sinking. Finally, Al rings for calibration 
purposes were fabricated in a similar way. Before each 
metal deposition, the developed resist was descummed in 
an oxygen plasma. The base pressure of our evaporator, 
which has never been used for magnetic materials, was 
below 5 X 10~^ Torr. On sample II, the first layer, an 
e-beam defined, 80 nm thick Au wire grid and bonding 
pads, was evaporated onto a Si substrate with native ox- 
ide from a source with unknown purity on top of a 1 nm 
Al wetting layer. A 50 nm thick AlO^^ film, patterned 
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using optical lithography and liftoff, was then deposited 
by atomic layer deposition (ALD). Rings and heat sinks 
similar to those on sample I were fabricated on top of the 
A\Ox [Fig. [T](b)], also without adhesion layer. 

Our dilution-refrigerator based microscope [13] em- 
ploys SQUIDs [1^ with an integrated 14 fim mean di- 
ameter field coil that is concentric with a 4.6 /im pickup 
loop [Fig. [T] (a)]. These loops can be brought to within 
about 1 /im of the sample surface. The field coil applies 
an ac field Ha (35 G amplitude at its axis for most of the 
data discussed here, corresponding to a field coil current 
/pc = 35 mA) to the sample, whose response couples a 
flux ^ SQUID into the pickup loop. A second, counter- 
wound pair of coils, located further from the sample, 
cancels the response of the SQUID to the applied field 
to within one part in 10^ [ij]. As the field coil cur- 
rent varies sinusoidally in time (with amplitude /pc ) , the 
SQUID response <1>squid is conveniently characterized in 
terms of its complex nth harmonics, <l>(^). We define 
(1)^1'^ +i(t)^^^ = ^^^V^FC and abbreviate ^1,2 = ^^^^ 
and (1)2 quantify the in-phase and out-of-phase linear re- 
spouse, ^ is proportional to the cubic component. 

Fig. [1] shows 2D susceptibility scans of both samples. 
For sample II, we took scans as shown in Fig. mc)-(e) 
at a range of temperatures and extracted the tempera- 
ture dependence [Fig. [2] (a), (b)] by averaging the indi- 
cated rectangular regions. For sample I, we averaged the 
complete /fc-^squid curves from many sinusoidal field 
sweeps at discrete positions as indicated in Figs. H^g), 
(j). This procedure is more sensitive and allows us to 
determine the difference between the full responses in- 
cluding nonlinearities [Fig. [2](d)-(g)] near and away from 
the metal. Figs. [Hh), (k) show 0i extracted from the 
response curves at each point. The T and frequency de- 
pendencies of the magnitude of the spatial variation along 
these line scans are included in Fig. [2](a)-(c). 

While the above measurements can only detect lateral 
variations of the sample response, the height dependence 
of the latter confirms that the observed signals reflects 
a response from the metal film. The height dependence 
over a third sample with a 1 /im thick SiO^c layer, grown 
at 1000° C with a wet process, and above Si with only a 
native oxide give evidence for a paramagnetic 1/T surface 
response that is about a factor 5 and 30 smaller than than 
that of the metal films, respectively [isj. The response 
from the ALD grown AlO^; film, which likely has a higher 
defect density than the thermal SiOa^, is comparable to 
that from the metal films [Fig. EJ^a)]. Comparison with 
the value of (j)2/4>i for superconducting rings [Fig. [2](c)], 
shows that the nonzero value of (j)2 is not an instrumental 
artifact. Nevertheless, the actual sample contribution to 
(j)2 is somewhat smaller than the raw values displayed in 
Fig. [lb) [H. 

The 1/T dependence and the paramagnetic sign of the 
susceptibility x(^) = Xi + ^X2 indicate that it originates 
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FIG. 1: (color online) (a) SQUID field coil and pickup loop, 
(b) Schematic of the layer structure of sample II: 1, Bare Si; 
2, Au with Al adhesion layer; 3, ALD deposited AlO^; 4, Au 
on AlOcc- (c)-(e) Sample II, linear in- and out-of phase signal 
(01, -02) and in-phase 3rd harmonic (0i^^) at 193 Hz, 43 mK. 
All numbers in this figure are in units of /i$o/mA and the 
response over bare Si has been defined as 0. (f) Scanning 
electron micrograph of a region of sample I, Au films on Si 
with native oxide, (g) Sample I: linear in-phase signal (01 ) at 
193 Hz, 27 mK of a region as shown in (f). (h) Line scans over 
the positions indicated in (g), at 25 mK and 111 Hz. Panels 
(i)-(k) are the same as (f)-(h), zoomed in on ring as indicated 
by the box in (g). The ring has a 2 /im diameter, 350 nm 
line width, and a connection to the 15 /im wide wire for heat 
sinking. The line scans in (k) were taken at 35 mK and the x- 
scan is offset for clarity. The temperature dependence shown 
in Fig. [2] (a) and (b) was obtained from line scans as in (h) 
and (k), or by averaging over the rectangles in (c). 



from localized spins. One can show that the 2:-field ema- 
nating from a film of thickness d with an isotropic linear 
response is B = jiiox^dtl/dz^ where H is the applied field 
reflected about the film (xy) plane. Using the measured 
pickup-loop-field-coil inductance and modeling H as the 
field of a thin loop leads to Xi,2<^ = 8 /im-mA/<I>o • 0i,2- 
For the films, this implies xi^ = 3 • 10~^ K to within a 
factor of 2. We estimate systematic errors of less than 
a factor 2 due to the simplicity of the model and uncer- 
tainties in the scan height. The response of the ring is 
consistent with this estimate within its somewhat larger 
calibration uncertainty. Comparing xi^ = 3 • 10~^ K 
to the susceptibility xi = Mo^(^Mb)^^(^ + 1)/3/cbT of 
dilute spins with number density n and total angular mo- 
mentum J leads to a concentration of 60 ppm for d = 100 
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FIG. 2: (color online) (a) Temperature dependence of the 
in-phase linear susceptibilities (01 ) obtained from scans as 
shown in Fig. [T] (see text) at 111 Hz or 193 Hz. Numbers in 
parenthesis indicate the sample. The downturn at the highest 
T can be attributed to the T independent diamagnetic bulk 
susceptibility of -3.4 • 10"^ for Gold O. (b) Out-of-phase 
component (02 ) from the same data sets, (c) Frequency de- 
pendence of 02/01 from sample I at 25 mK. The data from 
the superconducting ring characterize the phase shift due to 
the finite measurement bandwidth and show that the nonzero 
02 is not a measurement artifact, (d), (f) Nonlinear part of 
the response of the heatsunk Au ring (I), deposited directly 
on the Si substrate, at base temperature, (f), (g) Same for an 
isolated ring on sample H, deposited onto the AlO^ film. In 
(d), (f), only 01 (i.e. a line) was subtracted. In (e), (g), both 
01 and 02 (i.e. an ellipse) were subtracted. 



nm, g = 2 and J = 1/2, corresponding to an area density 
of 4 • 10^^ spins/m^. Because all our films have a similar 
thickness, we cannot distinguish whether the spin den- 
sity scales with the volume or surface area. For common 
magnetic ions with g^J{J + 1) ^ 35 [l^, the concen- 
tration would still be about 3 ppm, which is an order 
of magnitude larger than the specifications of our source 
material. For the materials investigated, the calculated 
equilibrium response from nuclear spins is several orders 
of magnitude smaller than the observed signals. 

The magnetic moment noise spectral density from a 
sample of volume V and susceptibility x that couples to 



a sensing loop as a dipole is S'm(^) = —2kBTx2{^)V/7rLJ. 
This form of the FDT was verified for a spin-glass, using 
a SQUID susceptometer similar to ours [l8|. In our case, 
the expected magnetic noise from the sample is much 
too small to be detected directly. However, Ref. H shows 
that a spin density similar to our estimates can explain 
the observed 1//" noise levels in SQUIDs. Integrating 
the dynamic susceptibility of fluctuating spins over the 
relaxation time distribution assumed in Ref. § leads to 
a value of 02 / that is also of the same order as our re- 
sults. Although we cannot prove that 0i and 02 are of the 
same origin, the similar T dependence, and consistency 
of the signs with a lag due to a finite relaxation rate, do 
suggest a direct connection. While it is also not a pri- 
ori clear to what extent our results apply to other metals 
including superconductors, the similar phenomenology of 
1// flux noise, our data from Au and Ag films, and recent 
measurements on superconducting devices also showing 
a 1/T susceptibility component \d\ indicate that all these 
effects are closely related. 

The frequency dependence of X2 and the nonlinearities 
imply a millisecond-scale spin relaxation time for some 
spins, which indicates weak coupling to the conduction 
electrons. On the other hand, we find that the linear sus- 
ceptibility of isolated metal rings saturates below approx- 
imately 150 mK [Fig. [2] (a)]. In such rings, the electrons 
are expected not to cool below that temperature because 
of heating by Josephson oscillations in the SQUID f\3\ . 
This observation indicates that the spins thermalize with 
the electrons rather than the lattice, which suggests an 
electronic relaxation mechanism 0] . 

One may thus suspect a connection with evidence for 
spin impurities in metallic nanostructures and at surfaces 
and interfaces obtained from transport measurements. 
Enhancement of superconductivity in nanowires in an 
applied field indicates pair breaking by spins [H, 
Weak localization measurements, which are a very sensi- 
tive probe for magnetic impurities [2l|, show that TiO^c 
adhesion layers for Au wires [22[ and native oxides on 
Cu films [23[ can cause spin-flip scattering. In contrast, 
our observation of similar susceptibilities from Au films 
with (wires sample I) and without Ti layer (rings I and 
II) shows that TiOx is not the dominant source of spins 
in our samples. From standard weak localization mea- 
surements for T > 300 mK on wires fabricated together 
with the samples discussed above, we find a dephas- 
ing rate l/r^ with a temperature dependence close to 
l/r^ oc T^/^, as expected for electron-electron interac- 
tion mediated dephasing The deviation from this 
power law behavior can be accounted for with 0.1 ppm 
of Mn impurities, if one allows the prefactor of the T^/^ 
term to be a factor 6 larger than theoretically expected. 
In typical weak localization measurements, the discrep- 
ancy between the theoretical and experimental prefactor 
is no larger than a factor two [2l|. However, if the un- 
usually large discrepancy in our case were due to spins. 
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their Kondo temperature would have to be larger than 
about 1 K in order to explain the observed increase of 
T(j) at low T. While we cannot rule out the existence of 
such spins, they cannot explain the susceptibility signal 
because their response would be quenched by the Kondo 
effect at low T. On the other hand, the measured > 1 
ns for T < 1 K sets a low upper bound on the spin flip 
scattering rate and thus exchange coupling and Kondo 
temperature of the spins contributing to the susceptibil- 
ity response 0. This upper bound does not necessar- 
ily rule out the RKKY coupling proposed in Ref. [Tlj . 
where it was suggested that the magnetic noise is due to 
RKKY-mediated spin diffusion. However, spin diffusion 
mediated by an isotropic interaction as considered in Ref. 
Uj conserves total angular momentum and thus the to- 
tal magnetic moment (assuming no ^-factor variations). 
Thus, the observation of a paramagnetic response from 
isolated rings, which are smaller than the pickup- loop 
[l2| and thus couple to our sensor mostly through their 
total magnetic moment, is inconsistent with this diffusion 
model. Nevertheless, anisotropic spin-spin interactions, 
such as dipolar coupling, could in principle determine the 
relaxation dynamics. 

We finally discuss a few anecdotal observations. We 
observed comparable values of (pi on two other samples, 
one similar to sample I, but on a Si substrate with an 
approximately 1 /im thick wet thermal oxide [lH, and 
one similar to sample II, but with Ag substituted for the 
top two Au layers. The Ag films showed a substantial 
spatial variation of the magnetic response. This inhomo- 
geneity could be due to an inhomogeneous surface oxi- 
dation or other chemical contamination. We find that 
the contributions from different layers on sample II are 
in general not additive, which most likely means that the 
spin population is concentrated at surfaces or interfaces. 
The nonlinear response and (j)2/4>i vary significantly be- 
tween different samples and different layers [Fig. Efd)- 
(g)]. Even though the nonlinearity seen in sample II is 
predominantly cubic, its magnitude and T dependence 
are inconsistent with the saturation of the equilibrium 
response at finite field. The relatively small ratios of 
the hysteretic and nonlinear components to (pi indicates 
that the majority of the spins contributing to the latter 
relax fast compared to the measurement frequency. The 
increase of ^2 with frequency is qualitatively consistent 
with varying slower than as observed in some 

SQUIDs 0, 
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In summary, we have measured the susceptibility of mi- 
cropatterned thin films. Different samples showed similar 
linear responses corresponding to an area density of un- 
paired spins on the order of 0.4 spins/nm^. The spins 
on our Au films appear to be weakly coupled to con- 
duction electrons. The out-of-phase component of the 
susceptibility gives direct experimental evidence for the 
hypothesis that the 1// flux noise seen in SQUIDs and 
superconducting qubits is due to fluctuating spins 
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